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ABSTRACT 

Super Star Clusters {M cc \ > 10 5 Mq) are the largest stellar nurseries in our local Universe, containing 
hundreds of thousands to millions of young stars within a few light years. Many of these systems are 
found in external galaxies, especially in pairs of interacting galaxies, and in some dwarf galaxies, but 
relatively few in disk galaxies like our own Milky Way. We show that a possible explanation for this 
difference is the presence of shear in normal spiral galaxies which impedes the formation of the very 
large and dense super star clusters but prefers the formation of loose OB associations possibly with 
a less massive cluster at the center. In contrast, in interacting galaxies and in dwarf galaxies, regions 
can collapse without having a large-scale sense of rotation. This lack of rotational support allows the 
giant clouds of gas and stars to concentrate into a single, dense and gravitationally bound system. 
Subject headings: ISM: clouds - open clusters and associations: general - galaxies: star clusters: 
general - galaxies: star formation 



1. INTRODUCTION 

Stars generally form in loose groups or embedded 
clusters in Giant Molecular Clouds (GMC), eac h cluster 
containing a dozen to many million of stars dTutukov 
19781: iZinnecker et all 119931; lLada k Ladall2003t iKroupa 
200a lAllen et all 120071) . These embedded clusters 
need not to be bound or radially well-defined stellar 
ensembles and the vast majority of them (~ 90% 
will disperse withi n about 10 Myr ([Tutukovl 1197 
lLada k Ladal I2003D . Only a small fraction of these 
clusters will hatch from the clouds as bound clusters, 
and after a consi derable loss in stellar numbers due 
to ga s expulsion (iGoodwinl 119971 iGoodwin k Bastiar] 



ent (IWeidner k Kroupall200i lFigerll200a IQev k Clarkd 
120051 : IKoeiJ 1200ft IWeidner et al.ll2010lT Surprisingiv. it 
is in the much smaller LMC that we find the largest 
stellar cluster in the local universe, 30 Doradus. At the 
center of the roughly 200 pc wide 30 Doradus Hll re- 
gion lies the "star-burst cluster" or super-star cluster 
(SSC) NGC2070/R136. It contains about half a mil- 
lion stars with an estimated total mass of ~ 2-10 5 Mq, 
extrapolated from the sa 50000 M Q in stars from the 
limited r ange of 2.1 and 25 Mq within an area of about 
150 pc 2 dSelman et all [l999l I Andersen et al.ll2009T ). On 
larger scales. iBosch et ~ 1 (120091 ) derive a dynamical 
mass of 4.5T0 5 Mq. It therefore qualifies as a possi- 
ble globular cluster precursor. The most-massive star- 
forming regions known in the M W have only s=s 50000 Mq 
( Wester lund 1. lBrandner1l2008fl to ps 75000 Mq (Arches 
and Cygnus OB2 iFiger et alj|2002HMartins et all 120081: 
Wolff et all 120071 : iNegueruela et alj|2008t IWeidner et al l 
20101 ). Though, as no strict definition of SSC's exist, 



1200ft IWeidner et al.ll2007t iBaumgardt k Kroupall2007l ). 
While these embedded clusters vary hugely in their 
number of stars, their physical sizes are all found 
to b e rather similar w it h R < 1 pc dTesti et al.l 
1991 [Gutermuth et all [20051: iTestor et all 1200a 



Rathborne et a l. 2006; Scheepmaker et al.ll2007f ). As the 



bulk of the galactic field star populations are probably 
made from dissolving embedded clusters dKroupal 
19951 lLada k Ladal 11991 l200l lAdams k Mversl 120011: 
Allen et aLll2007l ). studying the formation and evolution 
of the systems is very important. 

When comparing the star cluster populations of the 
Milky Way (MW) and the Large Magellanic Cloud 
(LMC), we see that although the stellar populations are 
similar, the LMC has much more massive clusters than 
have been found in the MW. The stellar initial mass func- 
tion (IMF ) do not differ beyond the expected statistical 
variation (Masscy 2003; Larson 2002: IWvse etldl 12001 : 
lKroupall2002t ). Furthermore, despite the lower metallic- 
ity of the LMC and other dwarf galaxies, the mass of 
the most massive stars does not seem to be any differ- 
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the Arches cluster clos e to the Galactic Cente r and NGC 
3603 (« 1.3-10 4 M , lHaravama et alJ 12008) are some- 
times called SSCs. But the most massive clusters in the 
MW appear to have masses 5 to 10 times smaller than 
30 Doradus. It is possible that comparable mass clusters 
do exist in the MW but are obscured. 

Massive or super-star clusters appear more frequent 
in interacting galaxies and based on our knowledge 
of the local group, in dwarfs galaxies. For example, 
in the interac t ing A ntennae galaxies (NGC 4038/4039) 
I Zhang fc Falll (|1999l ) found a large number of very mas- 
sive (> 10 6 M Q ) SSCs. In Tab. CD and Fig. Q] are shown 
the most-massive young (< 50 Myr) star cluster for a 
sample of local galaxies together with the star-formation 
rate (SFR) and the type of the galaxy. For the whole 
sample there seems to be no trend with galaxy type and 
most-massive star cluster but when restricted only to 
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TABLE 1 

Most-massive young (< 50 Myr) star clusters in galaxies 
of different types. 'SA' indicates spiral galaxies, 'SB' 
barred spirals and 'Irr' irregular and dwarfs galaxies. 



Irr/dwarfs 



Fig. 1. — Mass of the most-massive young (< 50 Myr) star cluster 
in nearby spiral, barred spiral and irregular (dwarf) galaxies (see 
Tab. [TJ. Filled symbols indicate mergers and galaxies with signs 
of recent interaction. 



non-interacting galaxies (open symbols) there is a clear 
trend with galaxy type. Irregular and dwarf galaxies, 
which are systems with less shear, tend to have more 
massive star clusters than non-interacting spiral galax- 
ies. But due to the small number of objects and large 
errors in the cluster masses, the statistical significance of 
the effect is low. 

This difference is even more striking when one consid- 
ers that dwarf galaxies have SFRs 5 to 10 times lower 
than do spirals, as can be seen by the mean values of the 
SFR and the maximum cluster mass for different types 
of galaxies in Tab. Q] A naive speculation would be that 
such galaxies should only produce lower-mass clusters, 
when in fact dwarfs appear to produce more massive clus- 
ters than (non-interacting) spirals. 

One main difference between dwarf galaxies, disk 
galaxies and interaction regions of galaxies is the amount 
of shear acting on GMCs in them. While disk galaxies 
are rotationally supported systems, resulting in relatively 
large shear forces on GMCs, dwarf galaxies show far less 
amounts of shear. Interacting galaxies with tidally driven 
structures are also likely to have low shear in the tidal 
arms and regions of interaction where SSCs are formed 
(see also § [2] on more details about shear in disk and 
dwarf galaxies). 

The presence of retrograde rotating GMCs (Blitz 1993) 
in spiral galaxies is sometimes seen as evidence that the 
formation and evolution of GMCs is not or only weakly 
connec t ed to the shear/rotation of spiral arms. In fact 
iDobbsl (|2008f ) has shown that spiral shock models of 
GMC formation produce both prograde and retrograde 
rotating GMCs due to the random nature of the coagula- 
tion process. Their internal shear is still directly related 
to that of the galaxy. 

In this paper we will examine the influence of different 
levels of shear on the fragmentation properties of GMCs 
with masses of 10 6 M Q . The model is described in § [2 
while § [3] shows the results which are discussed in § S) 



Name 


[10 4 Mo] 


SFR 

\Mq yr" 1 ] 


Type 


Ref. 


NGC 45 


3.6 


0.0254 


Irr 


(1) 


NGC 247 


8.9 


0.0364 


Irr 


(1) 


NGC 1156 


9.0 


0.186 


Irr 


(1) 


NGC 1569 


40.0 


0.037 


Irr 


(1) 


NGC 1705 


50.0 


0.03 


Irr 


(1) 


NGC 5585 


21.4 


0.0341 


Irr 


(1) 


NGC 7793 


16.5 


0.141 


Irr 


(1) 


IC 1613 


0.3 


0.00045 


Irr 


(2) 


LMC 


45.0 


0.2 


Irr 


(3) 


SMC 


36.0 


0.15 


Irr 


(4) 


mean 


23.0 


0.08 






NGC 300 


5.0 


0.0774 


SA 


(1) 


NGC 628 


20.4 


0.993 


SA 


(1) 


NGC 3184 


16.5 


0.395 


SA 


(1) 


NGC 5055 


11.9 


0.12 


SA 


(1) 


M 31 


15.0 


0.34 


SA 


(5) 


M 33 


4.0 


0.21 


SA 


(6) 


mean 


12.0 


0.35 






NGC 1313 


24.8 


0.426 


SB 


(1) 


NGC 2403 


9.4 


0.338 


SB 


(1) 


NGC 2835 


8.3 


0.0924 


SB 


(1) 


NGC 3521 


5.6 


2.19 


SB 


(1) 


NGC 7424 


28.4 


0.173 


SB 


(1) 


MW (Arches) 


7.7 


1.07 


SB 


(7) 


mean 


14.0 


0.71 






NGC 2997** 


104.0 


5.0 


SA 


(1) 


NGC 4038/4039* 


200.0 


8.3 


SA 


(8) 


NGC 5194* 


35.0 


4.75 


SA 


(1) 


Arp 220* 


900.0 


240 


SA 


(9) 


NGC 3621* 


38.6 


0.881 


SB 


(1) 


NGC 4258* 


30.8 


0.69 


SB 


(1) 


NGC 5236* 


38.6 


2.30 


SB 


(1) 


NGC 6744* 


10.7 


0.434 


SB 


(1) 


NGC 6946** 


85.1 


2.54 


SB 


(1) 


mean*** 


68.0 


3.1 







* Interacting or merging galaxy. ** Only i ndirect evidence for a 
recent interaction in N GC 2997 (Hess ct al. 2009|) and NGC 6946 
(Boomsma et al.l 12008). *** Arp 220 is not included in the mean 

values for intera cting galaxies. 

References: 1 : IWeidner et al l J2004T) ; ILarsenl ([2003. 12(3097); 2 : 
I Wvder et all 0200011: ILarsenl 1120021); 3: lA ndcrscn et al.l 12009T); 
Harris & Zaritsky (2009); 4: Harris & Zaritsky (2004); Sabbi ct al. 
I|2008|) ; 5: ICaldwell et all (120091 ); ITabatabaci & Bcrk huiisenl 
Mmf)l 6: IThilker et all IT2002T); iGrossi et all (12010) ; 7: 

I Weidner etall 1120101; I Robitaillc fc Whitney (2010)j 8: 

Izhang fc l^jTfll999l): Iwhitmore et all 112010); 9: I Wilson et al] 



(I2006T) ; 



2. THE MODEL 



The simulations of the gravitational collapse of GMCs 
under different initial conditions were carried out using 
a three dimensional (3D) smooth particle hydrodynam- 
ics (SPH) code, a Lagrangian hydrodynamics formalism 
(|Monaghanlll992l ). All clouds started from cold initial 
conditions in terms of both their thermal and turbu- 
lent energies being significantly subvirial and therefore 
collapse rapidly d ue to their s elf-gravity. A barytropic 
equation-of-state (|Larsonll2005[ ) is used with 7 = 0.75 for 
densities less than 5.5 x 10 -19 g cm -3 and 7 = 1.0 for 
densities above. The initial temperature is ~ 50 K and 
the minimum temperature reached is 7.5 K. The turbu- 
lent energies are w 0.06 that of the magnitude of the 
gravitational energies. The details of the initial condi- 
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tions are summarized in Table [5J In each case the gas is 
represented by 10 6 SPH-particles for the 10 6 Mq clouds. 
To m odel the star formation, sink particles (Ba te et"aLl 
fl9"95h are used which can grow through accretion of in- 
falling gas (SPH particles) and interact gravitationally 
with the rest of the simulation. The radius of the sink 
particles is r S i n k = 0.05 pc in all four cases. Each sink 
usually starts with about 50 to 100 Mq and can accrete 
up to ^5000 Mq over the length of the simulation. The 
sinks are therefore not to be seen as a single stars but 
rather as a small sub-clusters of stars. 

The calculations are evolved for about one free-fall 
time of the GMC. The free-fall time, tg, is the time- 
scale an object needs to collapse from its current radius 
into a single point u nder the influence of grav ity alone 
and can be given as (jBinnev fc Tremaineiri987D : 

1 / 3tt \ 1/2 

with G being Newton's gravitational constant and pgmc 
the mass density of the GMC. With a radius of 
-Rgmc = 50 pc and a mass Mqmc = 10 6 Mq the free-fall 
time is iff = 5.9 Myr in all cases considered here. 

The temporal evolutions of four model-GMCs are stud- 
ied, each with the same total gas mass (Mqmc = 10 6 
Mq) and radius (-Rgmc = 50 pc) but with four dif- 
ferent levels of shear. The first model has no shear, 
while the remaining three have shear levels correspond- 
ing to solid body rotation with angular velocities of = 
2 x 10" 15 , 5 x 10" 15 and 10" 14 rad s" 1 . The non-rotating 
model is taken to replicate the shear level expected in ei- 
ther a non-rotating dwarf galaxy or one in which the 
interaction with another galaxy produces regions of low 
shear in the tidally-induced spiral arms. For example, in 
the LMC a low shear value of Qlmc ~ 6 x 10~ 16 rad s _1 
is inferred from t he maximum of the ro tational velocity 
curve of its stars (jAlves fc Nelsonll2000D . 

Models 2-4 represent typical conditions expected in 
spiral galaxies. The pattern speed of the Milky Way's 
spiral arms is estim ated as Qmw ~ 10~ 15 rad s _1 
(jBissantz et al.l 120031 ) . In addition, spiral arms com- 
press the gas from different galactocentric radii which 
will increase the local shear rates by factors of 10 or 
more. For example, the velo city g radient in the Orion A 
GMC (M G mc ~ 2-10 5 M Q jBhtjlllMI) is about 0.3 x 
lO^rads- 1 (|Kutner et al.lll977tlBldiv et al.lll987[ ). For 
several other GMCs in the Milky Way very similar val- 
ues have been measured, e.g like Roset te, Mon Rl and 
W3 have . 3 to 0.6 x 10~ 14 rad s" 1 (jThronson et al.1 
119851: iBlitzl Il991h . Furthermore, values like the MW 
ones are found for GMC s in the local spiral galaxy M33 
(|Rosolowskv et al.l 120031) . All these GMCs have masses 
from several 10 4 to a few 10 5 Mq, roughly similar to the 
GMCs in our numerical study. 

The simulations are followed over ~6 Myr, where the 
free-fall time of the cloud is tg — 5.9 Myr (see eq. []}. 
Each model ran for about 2 month on the SUPA Altix 
computer of the University of St Andrews. 

The simulations presented here do not include any 
feedback from supernovae, radiation or ste llar winds. 
Thou gh, the effects o f ionizing radiation dDale et all 
I2005D and stellar winds pale fc Bonnellll2008D have been 
studied before. We note that the inclusion of these 



sources of feedback in the above models did not have a 
significant effect on the star formation rate or efficiency. 

3. RESULTS 

The evolution of the four models is shown in Fig. [2j 
The cloud evolves due to the internal turbulence, forming 
filamentary structures. These structures contain local 
regions which are gravitationally unstable and collapse 
to form sink-particles, localized regions of star formation. 
As the collapse continues, more sink-particles are formed 
and these join together in clusters which subsequently 
merge into one large super-cluster containing over 10 5 
Mq within 1 pc, and close to 4 x 10 5 Mq within 10 pc. 

At 4.3 Myr in Fig. [2] it can be seen that for increasing 
levels of shear the collapse and fragmentation produces a 
more distributed population rather than the highly con- 
centrated super-star cluster found in the no-shear run. 

All four models are evaluated at a time of about 
4.3 Myr (Fig. [2]) after the beginning of the calcula- 
tions. As the first sinks formed slightly after 1 Myr, 
some of the massive stars will have reached an age of 
3 Myr at which point they would explode as a supernova 
(|Mevnet fc Maederl l2003f ). We halt the calculations at 
this point as we do not include feedback in these models. 

As can be seen in Fig. [21 there is a strong dependence 
of the amount and central concentration of stars on the 
strength of the applied shear. At lower levels of shear, 
more sinks are formed and they are much more highly 
concentrated towards the center of the cloud. Figure [3] 
shows the temporal evolution of the central cluster core 
in each model. The mass density is calculated in a sphere 
of r corc = 1.0 pc around the center of mass of all sinks. 
The models diverge after 3 Myr as rotational support 
in the shear models halt the central collapse and core 
formation that occurs in the absence of shear. The no- 
shear model reaches average central densities in excess 
of 6 x 10 4 M Q pc -3 , up to 15 times larger than the runs 
with shear. 

The difference between the models with and without 
shear are more evident when considering the mass distri- 
bution of the resultant clusters and the distribution of lo- 
cal stellar densities (Figs.|4]and[5]). In the no shear runs, 
the cluster contains nearly 10 5 M© inside 0.1 pc whereas 
the shear runs have a much more distributed popula- 
tion without any significant central condensations. The 
total stellar masses in all cases, and hence the star forma- 
tion efficiencies are not drastically dissimilar with star- 
formation efficiencies (SFE — M sinks /MGMc) of 44%, 
38%, 26% and 18% for model 1 to 4, respectively. It 
is simply the distribution of the resultant stellar popula- 
tions which are very different. Also shown in Fig. [4] are 
five star clusters and OB associations in the MW and 
the LMC. Especially, the MW objects in the disk (NGC 
3603 and Cygnus OB2) seem to follow the high shear 
calculation quite well while NGC 2070 in the LMC fits 
the low shear one. The two other MW objects, Arches 
and Westerlund 1, are found in between the low and high 
shear run. Interestingly, neither of these objects are in 
the Galactic disk but Arches is very close to the Galactic 
center and Westerlund 1 is at the outer edge of the Galac- 
tic bar, both tidally very different regions compared to 
the disk. 

A stronger indication is the cumulative mass distribu- 
tion as a function of local stellar density measured as the 
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TABLE 2 

Initial conditions of the four models shown here. Mqmc is the initial mass of the GMC, J?gmc the initial radius of the 

GMC AND f! THE INITIAL ANGULAR ROTATIONAL VELOCITY, a AND DESCRIBE THE FOLLOWING ENERGY EVOLUTION PARAMETERS, a = 
^thermal /-^gravitational I fturb = ^turbulencc/-E'gravitational AND Plot = ^rotational /-^gravitational • ALSO SHOWN ARE THE MASS IN SINKS, 

A/4.3Myr> AND NUMBER OF SINKS, N i , 3Myr , AFTER 4.3 MYR. 



No. 


Mgmc 

Mq 


-Rgmc 
pc 


a 

rad s~ 


-l 


a 




Aurb 


/3rot 


-W4.3Myr 

Mq 


-W4.3Myr 


1 


lxlO 6 


50 







4x10' 


-3 


0.059 


0.0 


4.4xl0 5 


1451 


2 


lxlO 6 


50 


2x10" 


-IS 


4x10- 


-3 


0.059 


0.032 


3.8xl0 5 


1326 


3 


lxlO 6 


50 


5x10- 


-IS 


4x10- 


-3 


0.059 


0.158 


2.6xl0 5 


1100 


4 


lxlO 6 


50 


1x10" 


-14 


4x10" 


-3 


0.059 


0.574 


1.8xl0 5 


736 




Fig. 2. — Time-evolution series of four models. Each image shows a box of 60 times 60 pc. With increasing shear (model 1 to model 4), 
the amount and concentration of star- formation decreases. 
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time [Myr] p 10 [M /pc 3 ] 



Fig. 3. — Temporal evolution of mass density in a core with r CO re 
= 1.0 pc. The black dots mark model 1, the turquoise crosses 
model 2, the red boxes model 3 and the blue triangles model 4. 
After about 3 Myr of very similar evolution of all four models 
the calculation without shear collapses more quickly and to higher 
densities. It reaches densities which are about 3 to 15 times larger 
than those in the runs with shear. 
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Fig. 4. — The radial dependence of the accumulated mass for all 
four simulations after 4.3 Myr. The solid black line is for model 
1 (no shear), the long-dashed turquoise line for model 2 (some 
shear), the dotted red line for model 3 (intermediate shear) and 
the dashed blue line for model 4 (high shear). Also shown are 
several star clus ters in the MW and NGC 2070 in the LMC. The 
masses are from Wcidncr ct al. (2010). 



density of the ten nearest- neighbors of each sink shown in 
Figure [S] From this we see that approximately 10 5 Mq 
is in stars which have local densities from 10 6 to 10 s Mq 
pc~ 3 in the no shear run whereas in the highest shear run, 
less than 10 4 Mq is in stars which have local densities 



Fig. 5. — The dependence of the accumulated mass on the density 
derived from the ten nearest neighbors of each sink for all four 
simulations after 4.3 Myr. The solid black line is for model 1 (no 
shear), the long-dashed turquoise line for model 2 (some shear), the 
dotted red line for model 3 (intermediate shear) and the dashed 
blue line for model 4 (high shear). Clearly, the SPH-calculation 
without shear accumulates the most mass and this mass is very 
concentrated compared to the other three cases. 

approaching 10 5 Mq pc -3 . Clearly, the SPH-calculation 
without shear accumulates the most mass and this mass 
is very concentrated compared to the other three cases. 

For the analysis of the simulations, each time step 
the center of mass is derived by searching for the high- 
est mass density not of the individual sink but of the 
ten nearest neighbors. As can be seen in the right 
panel of Fig. [5] the local mass density reaches rather 
high values (> 10 8 Mq pc -3 ), which is to be expected 
when, e.g., twenty sinks of masses ~ 400 Mq are in- 
side 0.025 pc. Even in the relat ively small (~ 2000 Mq, 
iHillenbrand k, Hartmannl 119981 ) Orion Nebula Cluster 
(ONC), roughly comparable stellar densities are found. 
The projected radius of the 4 Trapezium stars in the 
ONC is about 10 arc-seconds, which translates to 0.02 
to 0.025 p c for assumed dista nces to the ONC of 400 
to 500 pc (jMuench et all 120081 ). The combined mass of 
these four systems (all four stars are binaries or higher- 
order multiples) is - 100 Mq (jPreibisch et all 119991 : 
iKraus et al.1 120091 ). Therefore, the stellar density is of 
the order of 10 6 Mq pc -3 . A second point is that below 
0.05 pc the gravity in the SPH calculation is smoothed 
and hence 2-body interactions do not occur and systems 
are stable where otherwise they need not to be. The sinks 
can cluster well within 0.05 pc and not feel each other 
directly and hence be stable. Therefore, while the local 
mass density might be overestimated at times in the cal- 
culations, the main point still holds that large amounts 
of mass are concentrated at high stellar densities in the 
no-shear run. 

4. DISCUSSION 

Our simulations show the formation of super-star clus- 
ters can depend strongly on the shear content in the pre- 



6 



collapse giant molecular cloud. Clouds which are only 
slowly rotating or not rotating at all, as to be expected 
in dwarf galaxies like the LMC and in interaction regions 
of colliding galaxies, can collapse monolithically into a 
single massive star cluster. Th i s is in agreement with 
the results by Es cala fc Larson! (|2008l ) who predict the 
precursors of globular clusters to form only in galaxies 
not stabilized by rotation. 

In contrast, in our calculations clouds in disk galaxies 
are more prone to fragmentation and form a system of 
smaller clusters or structures that could evolve to become 
OB associations or relatively large clusters susceptible 
to the tidal field of the galaxy with an extended halo 
of stars. Interestingl y, in HST observation s of massive 
extragalactic clusters [Maiz- A pcllaniz (20Ql|) finds a sim- 
ilar result. Four galaxies (NGC 2403, NGC 1569, NGC 
1705, LMC) of that study are in common with Tab. [TJ 
For these, the brightest, most compact clusters are all in 
dwarfs whereas the brightest objects in the normal spi- 
ral NGC 2403 are more extended, classified as so-called 
scaled OB associations. 

Although the results of the numerical calculations do 
not exclude the possibility of SSCs formation in disk 
galaxies, but they show that the presence of shear in disk 
galaxies acts to impede the formation of very massive 
clusters. It is therefore unlikely that the massive clusters 
which were the progenitors of the present day globular 
clusters have formed in the disk of galaxies. They might 
have formed either during an initial monolithic collapse 
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